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Abstract 

We investigate by means of Monte Carlo simulations the fully connected 
p-state Potts model for different system sizes in order to see how the static 
and dynamic properties of a finite model compare with the, exactly known, 
behavior of the system in the thermodynamic limit. Using p = 10 we are 
able to study the equilibrium dynamics for system sizes as large as = 
2560. We find that the static quantities, such as the energy, the entropy, 
the spin glass susceptibility as well as the distribution of the order parameter 
P{q) show very strong finite size effects. From P{q) we calculate the forth 
order cumulant g4{N,T) and the Guerra parameter G{N,T) and show that 
these quantities cannot be used to locate the static transition temperature 
for the system sizes investigated. Also the spin-autocorrelation function C{t) 
shows strong finite size effects in that it does not show a plateau even for 
temperatures around the dynamical critical temperature Tq. We show that 
the A^-and T— dependence of the a-relaxation time can be understood by 
means of a dynamical finite size scaling Ansatz. C{t) does not obey the time- 
temperature superposition principle for temperatures around To, but does 
so for significantly lower T. Finally we study the relaxation dynamics of the 
individual spins and show that their dependence on time depends strongly 
on the chosen spin, i.e. that the system is dynamically very heterogeneous, 
which explains the non-exponentiality of C{t). 

I. INTRODUCTION 

In recent years it has been recognized that the relaxation dynamics of supercooled liquids 
and the one of spin glass systems have many properties in common [^]-^. Of particular 
significance was the observation by Kirkpatrick, Thirumalai, and Wolynes P], that for certain 
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mean field spin glasses the equations of motion for the spin-autocorrelation function C{t) 
are formally the same as the ones that have been derived previously for the relaxation of 
particle density correlation functions of structural glasses p. That result thus opened the 
possibility to find one common description for these two classes of glassy systems: The 
spin glasses with quenched disorder and the structural glasses in which the disorder is not 
included explicitly in the Hamiltonian. 

One of the important classes of spin glasses models is the Potts glass |^J^-[TP|, which is 



a generalization of the Ising spin glass p0|-|2^, where the spins cxj can assume two values 



(cTj = ±1), to the case where each Potts "spin" can be in one of p discrete states, e 
{1,2, ... ,p}, p being an integer. Just as the Potts ferromagnet p8|-p0|] is a "workhorse" 
for the statistical mechanics of phase transitions, since it has served to test many methods 
and to exemplify many concepts about the subject, one can expect that the Potts glass will 
play a similar role for the study of the glass transition in systems with quenched disorder 
or the properties of liquids close to their glass transition temperature, since the possibility 
to change p allows to describe different glass transition scenarios, e.g. from a continuous 
transition to a discontinuous one. 

It is important to note that in the Potts glass it is well established that for p > 4 one 
has both a dynamical transition at a temperature T^, where the relaxation time associated 
with C{t) diverges, and a static transition at a temperature Tq < Td, where a glass order 
parameter appears discontinuously, accompanied by a kink in the entropy (as well as in 
the internal energy). In contrast to this knowledge, for the structural glass transition the 
existence of an underlying static transition, although proposed long time ago pl] , ^ , is still 



an open question p3|-p^ . Thus a study of Potts glasses should help us to understand also 
better structural glasses and hence is a welcome addition to the large efforts made to identify 
which structural features distinguish the solid glass from the liquid from which it was formed 
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Besides this interest in the Potts glass as a possible prototype model for the structural 
glass transition, it can also be considered as a coarse-grained model for orientational glasses 
42-^. Experimentally these systems are created by random dilution of molecular crystals. 



which has the effect that at low temperatures the quadrupole moments of the molecules 



freeze in random orientations If the crystal anisotropy singles out p discrete preferred 



orientations (e.g. the 4 diagonal directions in a cubic crystal), a Potts glass model with 
p states may give a qualitatively correct description of the system. And, last but not 
least, the Potts glass model of course completes our knowledge about the different types of 



phase transitions and ordered phases that spin glasses can have p2|-p^, which provides an 
additional motivation for the large activity in this field |^,0-|T^, 

In this work, we use Monte Carlo simulations to study the Potts glass model. Our goal is 
to clarify to what extent this interesting and nontrivial mean field behavior which is known 
exactly in the thermodynamic limit, N ■—>■ oo, can be seen for finite N. In addition, we want 
to elucidate the dynamical behavior of the model in greater detail than has been done so 
far and thus help to clarify the reasons for non-Debye relaxation in glassy systems. 

In the present paper, we shall first define the model and introduce the quantities that will 
be investigated (Sec. 2). In Sec. 3 we summarize what is known about the static behavior 
of the model and describes our pertinent numerical results. Sec. 4 is then devoted to the 
dynamical properties in the high temperature phase, i.e. above To, and the finite size 
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behavior at T^, while Sec. 5 is concerned with the dynamical behavior of small systems 
in the low temperature phase. By following the relaxation of individual spins we are able 
to study also the "dynamical heterogeneities" ^51 , ^6 
summarizes our conclusions. 



at low temperatures. Sec. 6 finally 



II. MODEL AND SIMULATION METHODS 



In this section we define the Hamiltonian and the observables that we consider in this work. 
Subsequently we give the details on the simulations. 

The Hamiltonian of the p-state mean-field Potts glass of interacting Potts spins ai 
{i = 1,2, . . . , N) that can take p discrete values G {1,2, ... ,p} is defined as 



(1) 



The "exchange constants" (bonds) Jij are quenched random variables which we assume to 
be distributed according to a Gaussian distribution P{Jij) 



■ exp 



2{Ajy 



V2n{AJ) 

The first two moments Jq and A J are chosen as follows: 

Jo = [Jijl. = Jo/{N-l), Jo = 3-p, 



(2) 



(3) 



(^^)'-[4]..-i^dL = V(iv-i), (4) 

where [. . .]av denotes an average over all realizations of disordered bonds (while thermal 
equilibrium averages will be denoted as (...)). The scaling of the parameters Jo, A J with 
N was chosen such as to ensure a sensible thermodynamic limit both for the spin glass 
transition and for the ferromagnetic transition (which occurs for a certain range of values of 
Jo). Note that for p = 2 Eqs. (|T]) - (^) simply reduce to the Sherrington-Kirkpatrick (SK) 
model of a spin glass [^. We also mention that the term ^ J^j in Eq. (P is only included 

i<j 

for convenience, since it makes the mean energy of each system, i.e. for each realization of 
the disorder, go to zero for T oo. For the present choice of the parameter AJ, the spin 
glass transition in the replica-symmetric approximation, within which one finds a second- 
order transition for p < 6, occurs at a temperature = 1 |]^. (Here and in the following 
we set Boltzmann's constant ks = 1). While for some choices of Jq the system exhibits a 
transition to a standard ferromagnetic phase at a temperature Tp > T^, for our choice of 



parameters Tp falls far below and hence ferromagnetic order is of no of interest here ||7|JTl 
There exists in addition a second transition to a different type of spin glass phase (some- 
times called "randomly canted ferromagnetic phase") [0,^, at a transition temperature T2 
which is given by 

T2 = (p/2- !)/(!- Jo). (5) 
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For the choice of Jq given in Eq. one thus finds T2 = 1/2. Hence this choice ensures 
that at the temperatures of interest in the present study, T > 0.7, any effect of this second 
transition on physical observables should be negligible. 

For defining observables like the magnetization, the glass order parameter, time- 
dependent spin autocorrelation functions, etc., it is useful to choose a representation for 
the spins that takes into account the symmetry between their p possible states. This can be 



achieved by the so-called "simplex representation" [29|,0 in which the p states correspond 
to {p — l)-dimensional vectors Sx pointing towards the corners of a p-simplex, i.e. 



•S'a ■ Sx' = {p6xx' - 1) with A, A' = 1, 



(6) 



In our study we consider static as well as dynamical observables. Static quantities include 
the energy per spin. 



e = [{n)U/N, 



(7) 



the spin glass susceptibility Xsg, and the spin glass order parameter distribution function 
P(g). For defining a spin glass order parameter, we follow the standard method used for 



Potts glasses |jT2|JT^,|T9| to consider two replicas a and /? of the system, i.e. two systems that 
have identical bond configurations, and to make for each of them an independent Monte 
Carlo simulation. The order parameter tensor is then defined as 



1 ^ 

i=l 



1,2, 



(8) 



In an equilibrium simulation of a finite system with no external fields that couple to odd 
moments of the order parameter the symmetry is not broken. Hence it is useful to consider 
the root mean square order parameter q defined as p|, p!2|Jl6| , p!9[ 



p 



1 



fJ,,U=l 



1/2 



(9) 



and by calculating a histogram of q, i.e. by taking first the thermal average and then the 
average over the disorder, one can estimate the above mentioned distribution P{q). The 
second moment of this distribution is related to the reduced spin glass susceptibility Xsg- 



N 



Xsg 



p 



N 



p — 1 



(fP{q)dq. 



(10) 



(Below we will discuss the relation between xsg and the standard spin glass susceptibility 
XsG^ see Eqs. (^TJ) and (p6D.) If there is a second order transition to a spin glass phase, xsG 
should show a divergence at the critical temperature. Further interesting quantities related 
to the distribution P{q) are the reduced fourth-order cumulant |T2| , p!6| , p!9 



{p-l) 



1 + 



_j ML 



(11) 
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and a quantity called the Guerra parameter 



sav 



The reason for introducing these ratios of moment is that they are useful in the context 
of finite size scaling analyses of phase transitions. They are defined such that for N = oo 
they are zero in the disordered phase and nonzero in the ordered phase. In the finite size 
scaling limit the curves gA{T), or G{T), for different system sizes should intersect in a 
common point at the static phase transition point. In particular, G is a measure for the lack 
of self-averaging. 

To study the dynamical properties of the system, we will mainly focus on the autocor- 
relation function of the Potts spins, 

1 ^ 

Cit) = j^^^^ ^[{S^{t') ■ S,{t' + t))U . (13) 

Note that in thermal equilibrium C{t) depends only on the time difference t, i.e. it is 
independent of the second argument t' occurring on the right hand side of Eq. (p!3D . In 
practice, for the Monte Carlo sampling using the Metropolis algorithm [Q, the thermal 
averaging (■ ■ ■) is a time averaging over the initial times t'. 

We have also considered a rotationally invariant order parameter time-displaced correla- 
tion function CRiif) which is defined as 



with 



m 



Cmit) 



1 V^' 1 ^ 

— J2 ^nt) f and r{t) = ^ $^(S^,)'^(t)(5,)'^(0). (15) 

u,v=l \ i=l 



mi 



(14) 



Note that q^'^ is not the same quantity as q'^'^ defined in Eq. (|^), since the latter involves 
two replicas a and (3. However, for t — oo the thermodynamic averages of the two quantities 
are the same. This means that in this limit also q and q, from Eq. (P), are the same. 
Furthermore we mention that the expectation value {q{0)) occurring in Eq. ([T^ is equal to 
1, as long as there is no ferromagnetic ordering of the system. It is important to realize that 
CRi{t oo) is not zero, if is finite. This follows immediately from Eq. (0), since for 
large t the quantity q^'^{t) is of order and g(t), as finite sum over such quantities, is 

hence positive and also of order 1/^/N. From Eq. (|TD|) one also concludes that Cjij{t oo) 
is of order ^JxsG/N. 

In our simulations, we have investigated 5 different system sizes, = 160, 320, 640, 1280, 
and 2560. The number of samples used to approximate the quenched average [. . .]^^ over 
the bond disorder was 500 for = 160, 200 for A^ = 320, 100 for A^ = 640 and 1280 and 
between 20 and 50 for A^ = 2560 (depending on temperature). At not too low temperatures, 
T > 1, the straightforward Metropolis algorithm was implemented [Q, picking spins at 
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random and choosing randomly an orientation for them as a trial configuration. Depending 
on the energy difference AE between the old configuration and the trial configuration, the 
trial configuration was always accepted if /S.E < 0, else it was accepted with probability 
P = exp{—AE/T). For temperatures T < 1 the equilibrium configurations were generated 
with the "parallel tempering" technique P9|-pl|] . These equilibrium configurations can both 



be used to study the static properties of the model and as starting configurations to study 
the usual Metropolis dynamics in equilibrium and thus to calculate C{t), although at very 
low temperatures the relaxation is so slow that one cannot follow the complete decay of C{t) 
to zero. The total computing time (in units of single Pentium II processor with 400MHz) 
used for this study was of the order of 10 years. 

III. STATIC PROPERTIES OF THE 10-STATE POTTS GLASS MODEL 

In this section we first discuss the analytic results for the static properties of the model in 
the thermodynamic limit. Subsequently we compare them with the results of the simulations 
for finite N. 

If one calculates the free energy of the model given by Eqs. (|l]) - with the replica 
method |^D|-^, (without allowing for replica symmetry breaking), one obtains, depending 
on the value of Jq but independent of p, either a transition to a spin glass phase (where 
the spin glass order parameter go is nonzero) or to a ferromagnet (with a spontaneous 
magnetization mo). Within this approach and close to a critical temperature Tg the free 
energy density f{qo,mo) can be written as follows 0, 

- /(go, mo)/T = 2^^' ( ^ ~ 7^ ) ^0 + g^^o + ^'^mml + -UmiriQ + -u'qQml + ... (16) 

where with our choice of units Tg = 1, and r', u, r^, Um, u' are constants that are of no 
interest to us here. If the parameters are chosen such that the transition that occurs at Ts is 
to a spin glass phase, it is found that the order parameter distribution P{q) is a ^-function 
whose position depends on T (we consider here only the case m > 0): 

P{q)=6{q), for T>Ts (17) 

P{q)=Siq-qo), with go = ^ - , for T < T,. (18) 

Note that the first term on the right hand side of Eq. (|16]) can also be interpreted in 
terms of the spin glass susceptibility xsg as \Xsg1o Q Thus close to Tg one finds for xsg a 
Curie- Weiss law: 

XsG=[2r'{l-TjT)]-\ T>Tg. (19) 



^For the definition of a proper conjugate field to define the spin glass susceptibility, see references 
|22|,BT| 
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The coefficient r' in Eq. ([I6D is given by [0 (a result coming from the expansion in q of the 
free energy close to Tg) 



P 



1 /T. 



1 



T 



(20) 



2 \T 

We thus find for the susceptibility xsg around Tg (remember T, = 1 in our normalization): 



XsG 



2(p 



T 



(21) 



The difference between the standard spin glass susceptibility and the reduced one defined in 
Eq. ([T0| ) is just the factor (p — 1), related to the susceptibility of a system of non-interacting 
spins. 

It is well known that, if one allows for replica symmetry breaking, the prediction 

that there is a second order transition to a spin glass phase remains valid only if p < 4. For 
p > 4, a new type of first-order transition to a glass phase is predicted to occur at a temper- 
ature To which is higher than T^. Although at Tq the order parameter jumps discontinuously 
from zero to a value go > there is no latent heat involved in this transition. Instead of 
Eq. (|18D the order parameter distribution acquires now a double-5 function structure [^|9|,pA| , 



P{q) = [1 - w{T)]S{q) + w{T)5{q - q,{T)), T < T,. 



(22) 



with w(T) = 1 — T /Tq for T ^ Tq . While it is possible to calculate go and Tq analytically 
for p ^ 4+, see e.g. |]Tl| , 



(23) 



go = -(p - 4) + o{p - A)\ To~TsCx{p- 4)^ + o{p - 4f 



for integer p > 4 the correct values for go and Tq can be obtained only numerically [0. E.g., 
for our case of p = 10 the predicted values are 



To = 1.1312 and go(To 



0.452 



(24) 



In the disordered phase, the internal energy per spin e and entropy per spin s are given by 

EM 



P-IT, 
"2 T' 



Inp — 



P-1 fTs^ 
T 



(25) 



and a high temperature expansion gives |24 



T>To. 



(26) 



Within the replica symmetric Ansatz these expressions are correct for T > Tg. If one 
allows for replica symmetry breaking they hold only for T > Tq. Although no explicit 
analytic expression for e(T) and s(T) are known for T < Tq their value can be calculated 
numerically |l^ . 
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Finally we mention that within the replica Ansatz (symmetric or broken) neither Ts nor 
To depend on the choice of Jq from Eq. (|^), provided that the transition temperature Tp 
from the disordered phase to the coUinear ferromagnetic phase, discussed in section |I|, is 
not above Tq [^|11|. Tp is given by the following equation, for arbitrary p 



1 + Vl + 2(p-2)/j2. 



(27) 



DeSantis et al. used Jq = ^(4 — p) in which case Tp = = 1, independent of p. 
However, this case is rather special since then the transition temperature T2 to the randomly 
canted ferromagnetic phase, discussed in the Introduction, coincides with Tg = 1, as can be 
seen from Eq. (^ . For our choice of Jq = 3 — p we have instead T2 = 1/2 for all p and 
Tp = 8/(7 — \/65) ~ 0.531 for p = 10. Thus with this choice we hence make sure that the 
ferromagnetic fluctuations are still very small at T^ = 1, even if the system size is rather 
small. 

We now present our numerical results and compare them to the analytical predictions 
that we just discussed. Fig. |l|a shows a plot of the internal energy versus inverse temperature, 
over a wide temperature range (0.7 < T < 2), but still clearly above the temperatures Tp 
and T2 (remember that l3p = l/Tp ^ 1.88, P2 = I/T2 = 2.0). Also included are the 
theoretical predictions for ^ 00 obtained within the replica symmetric and one-step 
replica symmetry breaking theory, respectively. This figure reveals unexpectedly large finite 
size effects over a broad temperature regime in that, e.g., the energy for A^ = 640 coincides 
with the asymptotic result only if /3 ^ 0.6. For (3 > 0.6 clear deviations from the asymptotic 
solution are visible, which are larger for smaller A^. The numerical data for finite A^, in 
the range accessible to our work, reveal only a smooth crossover from the regime of the 
disordered high temperature phase to the regime of the low temperature glass-like phase, 
and no indication of the kink at Pq, predicted by the one- step replica symmetry breaking 
theory, is yet visible. As expected even for N ^ 00, there is no effect of the dynamical 
transition at Pd = ^/Td on static quantities like the energy. 

Of course it is also of interest to study how at fixed temperature the energy eAr(T) 
converges to its asymptotic limit eoo(T). Fig. |l]b shows that the energy difference eAr(T) — 

(T) scales like A^~^ both for temperatures above the static transition temperature Tq and 
for temperatures below Tq, while at Tq a different law, (oc A^~^/^), seems to hold (see inset). 
(Note that we plot here data for To instead of Tq since we have simulated more system sizes 
at this temperature. However, since the two temperatures are so close to each other this 
difference should not matter for the system sizes investigated here.) 

We have also calculated the temperature dependence of the entropy s(T). This was done 
by a thermodynamic integration of the free energy /: 



s{P) = Pe{P)-Pf{P)+ r e{P)dp . (28) 

We have used (3 = 0.5, a temperature at which our data is no longer sensitive to finite 
size effects and thus the replica solution is valid, so that we can use for the free energy 
(3f{(3) = (3e{(3) - s{(3) the mean-field value -9/16 - log(lO) (see Eq. (|5D). The integral 
over e was done by using a spline interpolation of our data, with 180 points for A^ =160, 
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320, 640 and 100 points for =1280, 2560. We think that alternative methods to calculate 
e(T), such as re-weighting techniques or methods to directly sample the density of states 
|[48|| , would not bring a significant advantage in our case. The results are shown in Fig. |I]c 
from which we recognize that s(T) shows similar finite size effects as the energy e(T). 
From Eq. (|25|) we see that the entropy at the static transition is 

s(To) = In 10 - ^{Ts/ToY ^ 0.544, i.e. s(To)/s(T = oo) ^ 0.236. (29) 

Thus we see that, while at the static transition temperature Tq the entropy has decreased 
to less than a quarter of its high temperature value, it is clearly nonzero (and nonnegative, 
of course). In supercooled liquids one often extrapolates the temperature dependence of 
the entropy (minus the vibrational entropy of the crystal) to zero and uses this to define 



the Kauzmann temperature Tk If one proceeds in the same way with the current 

model to obtain a "Kauzmann temperature" Tk where the entropy of the metastable high 
temperature phase vanishes, one obtains from Eq. 



Tk/T, = ^(In 10)-^/2 ^ 0.9885 (30) 

which is even below the "true" metastability limit = 1 of the disordered phase, where 
the (extrapolated) static glass susceptibility is divergent. (Note that the proximity of Tk 
and Ts happens accidentally for p = 10. E.g. for p = 5 the general result 0| Tk/Ts = 

1 /2 

[i(p— l)/lnp] implies Tk/Ts ~ 0.7882.) Also a strictly linear extrapolation of s{T) 
from Ts or Tq down to a temperature where this extrapolation would vanish does not give a 
meaningful result, of course. 

These results show that the idea to locate the static glass transition temperature by an 
extrapolation of the (configurational) entropy function s{T) in the disordered high temper- 
ature phase to s(T = Tk) = [^,0 can be completely misleading, even for a mean- field 
model that does indeed exhibit both a dynamical transition (at Td) and a static transition 
(at To). While Tk is always lower than Tq, it does not coincide with the stability limit of the 
metastable high temperature phase, and thus lacks any physical significance. Since for the 
case of polymers theories were formulated that suggest that Tk is the static glass transition 
temperature it is interesting to note that a simulation study of the glass transition in 
the framework of the bond fluctuation model found a decrease of s from its high temperature 
value to about 1/4 of this value, when T is lowered, but that subsequently the curve s(T) 



vs. T bends over and a well-defined does not exist |54]. The "configurational entropy" 



estimated by Gibbs and DiMarzio ||32| was simply the total entropy of their lattice model of 
polymers, just as the total entropy of our model shown in Fig. |l]b. Thus this demonstrates 
that the calculation of Tk in this way is most likely wrong. Of course, this "naive" way to 
define the Kauzmann temperature via the vanishing of the (configurational) entropy s(T) 
should not be confused with the approach of defining a "complexity" P,P^|-p7|,|5^] . In that 
approach one determines the number of basins in the free energy and defines Tk as that 
temperature at which this number starts to become exponentially large. 

Fig. P shows the reduced spin glass susceptibility as a function of the squared inverse 
temperature. This representation is adapted to the theoretical temperature dependence of 
this quantity, see Eq. (EB), which predicts at T~^-law at high temperatures. As expected 
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already from the behavior of the internal energy, even far above T, and Tq the convergence 
to the thermodynamic limit is rather slow. Unfortunately the analysis of the finite size 
behavior of xsG is not straightforward, as we will show in the following. For T < Tq and 
N ^ oo Eqs. (0) and (H) imply 

— 1 _ P - 1 1 fo-i\ 

iVan^ (l-T/nV ^'^ ^ 



since for T < Tq our definition for Xsg, Eq. (|TDp, simply picks up a contribution due to the 
nonzero spin glass order parameter go- As a result, Xsg for — oo should follow the straight 
dashed line in Fig. that represents the replica-symmetric solution for all < Tq~^, while 
for T^^ > Tq"^, Xsg simply converges towards the abscissa. This singular behavior of x^q 
is explained further in the inset, where we have added to Xsg term (Ts/T)'^. This sum 
gives unity for T > Tq and {Tg/Ty for T < Tq, as can be seen from Eq. (pG]). For very 
large but finite A^, Xsg ^ < exhibits a Curie- Weiss type divergence at Tq, but the 
amplitude of this effect is only of order 1/N, as can be seen from Eq. (|3T|). In order to 
analyze the finite size rounding of this singularity, one must consider that for finite the 
5-functions in Eq. ( ^2] ) are broadened into peaks of finite height and nonzero width 0. Our 
simulation results for P{q), see Fig. |^, do indeed give evidence that a second peak at go 7^ 
develops, distinct from the peak at small g that exists also in the high temperature phase. 
However, the statistical accuracy of -P(g) is not very high due to the well known fact that in 
the ordered phase this quantity is not self- averaging P^|-p^], and the number of realizations 



of the random couplings that we were able to study is insufficient to overcome this problem. 
Hence we are currently not able to do a proper analysis of the finite size effects of -P(g), see 
Fig. and thus cannot make a finite size analysis of Xsg- 

From Fig ^ we see that even in the high temperature phase, i.e. T > Tq, we see a peak 
in P{q) at a finite value of g. That this is, however, a finite size effect is demonstrated in 
Figs. and ^ where we plot -P(g) for different system sizes and show the first moment of 
P(g) as a function of N, respectively. We see that for temperatures well above Tq the first 
moment vanishes like iV~^/^. However, close to Tq this type of extrapolation would give a 
finite value of the moment. If instead an extrapolation with A^~^/^ is done, see inset, one 
finds again as expected that the moment vanishes. Note that for the second moment of P{q) 
we would have at T = Tq again a scaling of the type N~'^^^, as we found for the case of 



the energy. It is interesting to note that Parisi et al. ||70|], making use of replica symmetry 
breaking scheme, were able to calculate the finite size scaling exponents for the Ising spin 
glass, and obtained that [(g*')] scales like N~''^^ at the critical temperature (consequently 
one has that e scales like A^~^/'^). It is thus possible that the same kind of scaling holds also 
for the Potts glass, although one has to keep in mind that both the nature of the transition 
and of the replica symmetry breaking pattern is different. 

While the results shown so far demonstrate a rather encouraging qualitative consistency 
between the theoretical predictions and the numerical data, our results for the fourth order 



phenomenological attempt to describe the finite size behavior for the glass transition of Potts 
models has been made in [16|, but this approach is not followed up here, since it is not consistent 
with Eq. (22) in the limit to ^ oo. 
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cumulant g^, Eq. (pTD, and the Guerra parameter G, Eq. ([121), are clearly rather worrisome 
(Fig. It is seen that the three curves for (74 (iV,T) have a rather well-defined common 
intersection point at T ?a 1.31, and the three curves for G{N,T) have a rather well-defined 
common intersection point at T ?s 1.24. Motivated by the standard knowledge and experi- 



ence with finite size scaling at second order |^J5^ and first order phase transitions, 

such intersection points are commonly taken as estimates of the transition temperature 
2^^7[, However, comparison with the exact result for Tq, Eq. (P^, shows that these inter- 



section points are spurious, and hence cannot be taken as accurate estimates of Tq. This 
already is obvious from the data alone, because the two quoted temperatures are not in 
mutual agreement. We think that in reality neither the three curves for g4{N,T) nor those 
for G{N, T) intersect at a unique temperature. Given the relatively large error bars of the 
data, they only can define temperature intervals ATg4, ATg, in which the three intersec- 
tion points fall. As N ^ 00, presumably all temperatures of these intersections converge 
(slowly!) towards Tq. Since Tq falls distinctly outside the above intervals, this method of 
searching for intersection points, which is so successful for locating phase transitions in pure 
systems P8| , |55| , |56| , is a complete failure here. We emphasize this problem so strongly, be- 
cause such techniques are commonly used for studying phase transitions in systems with 
quenched disorder [^. Again we stress that an analytical guidance for the description of 
the finite size rounding of first order glass transitions would be very useful. 



IV. DYNAMICAL PROPERTIES IN THE HIGH TEMPERATURE PHASE 

In this section we briefiy review the theoretical predictions for the relaxation dynamics of 
the spins. Subsequently we compare these predictions with the results from the simulations. 

The theoretical results of Kirkpatrick et al. show that the Potts glass with j9 > 4 states 
has a "dynamical transition" at a temperature To > Tq, where non-ergodicity sets in 0. 
For T <Td, the spin correlation function C(t), defined in Eq. ([l^) , no longer decays to zero 



but gets stuck for t — > 00 at a nonzero value qEA{T), with |]T3 



Td = I.U2, QEAiT = Td) = 0.328 . (32) 

The details of this transition from ergodic (for T > where C{t ^ 00) = 0) to non- 
ergodic behavior (for T < Tjj), as well as the time dependence of C{t) for temperatures 
around are in fact described by equations formally analogous to equations proposed 
for the structural glass transition by idealized mode-coupling theory 0. The qualitative 
behavior of various quantities expected for ^ 00 is sketched in Fig. ^. Note that for 
T > Td and T < To we have go = Qea- In the temperature range Tq < T < Tjj we have, 
however go = and qEA > 0. 

In Fig. ^ we show the time dependence for the spin auto-correlation function C{t) for 
N = 1280 and all temperatures investigated. Here and in the following we will measure 
time in units of Monte Carlo Steps (MCS), i.e. the average number of updates per spin. 
Surprisingly we see that even for this rather large system size there is not yet any clear 
evidence for the development of a plateau for temperatures around T^. Note that in the 
thermodynamic limit this function should at T = decay to Qea, i-e. the horizontal line. In 
contrast to this our system with a finite size is always ergodic, since the free energy barriers 
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separating the various "valleys" in phase space remain finite at all nonzero temperature. Of 
course every finite system with no hard-core interactions is in principle ergodic. However, 
e.g. in structural glasses it is found that even a few hundred particles are sufficient to show 
a pronounced (effective) ergodic to non-ergodic transition. Thus it is rather astonishing that 
for the present model the finite size effects are so strong that even for = 1280 and at 
T = Td the existence of a plateau can hardly be seen. 

Also the time dependence of CRi{t), Eq. (0), shows strong finite size effects, as can be 
seen from Fig. ^d. We see that, contrary to naive expectation, the long time limit of CRi{t) 
is not zero but a finite constant. This constant depends on temperature and below we will 
discuss its origin and its dependence on system size in more detail. At any rate, we see that 
also this correlation function does not show a plateau even if T is close to To and hence we 
conclude that also C/j/(t) converges only very slowly to its behavior in the thermodynamic 
limit. 

In order to discuss the system size dependence of the correlation functions in more detail 
we show in Fig. |^ C{t) and Cjij{t) for different systems sizes and two temperatures. From 
Fig. 1^ we see that, at high temperatures, C{t) shows basically no system size dependence. 
For low T, however, the relaxation becomes quickly slower with increasing system size and 
also the shape of the curve changes noticeably. But even at the largest system sizes accessible 
at this temperature we are not able to see a clear two step relaxation as one would expect 
for a sufficiently large but finite system. 

The A^-dependence is different in the case of Cmit), Fig. |^. Here we see that even at 
high temperatures the correlation function depends on the system size. This is in agreement 
with the arguments given in the context of Eq. (JT^) that Cmit oo) should scale like 
l/\fN . That one actually find this size dependence is shown in the inset of Fig. 1^3. Instead 
of studying the function CRiit) one could of course try to consider the reduced normalized 
function (pit) = [CRj(t) — Cmit oo)]/[C/?/(0) — Cmit oo]. However, also this type of 
correlation function has its problems since on one hand the final asymptote CRi{t cxo) is 
only known to within a certain statistical error, and on the other hand it shows finite size 
effects at high temperatures at shorttim.es, i.e. where CRi{t) is independent of A^. In view 
of these problems with CRi{t) we will in the following focus on C{t) only. However, this is 
not a serious restriction, since in the thermodynamic limit these two functions should show 
at low temperatures the same time dependence anyway. That this is indeed the case for the 
simulations can be inferred from Fig. ^ where we show a parametric plot of CRi{t) versus 
C{t) at Td. We see that with increasing system size the curves do approach the diagonal, 
as expected. 

We now address the temperature and A^ dependence of the relaxation time r of the 
system. One possibility to define r is 

C{t = t) = 0.2. (33) 

Note that although the value 0.2 is somewhat arbitrary, it is a reasonable choice. The 
only important thing is that it is significantly less than the height of the plateau in the 
thermodynamic limit, qsAiT = Td), cf. Eq. (P^). (If we would define a time r' as C{t = 
t') = 0.5, on the other hand, r' would be finite also below Td, and even below Tq, until 
(1ea(T) has increased up to qea = 0.5, due to the temperature dependence of the order 
parameter, see Fig. ^ja.) 
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Since for N oo the dynamics of the model should be described by (idealized) mode- 
coupling theory |Q, we expect that t(T) shows a power law divergence at P], 

r oc {T/Td - 1)"^, ^ oo , (34) 

where A is an exponent which is non-universal (i.e. model dependent), but typically not 
very different from A ^ 2. In order to test the validity of Eq. (Bl), one can plot r-i/^ for 
a reasonable trial value of A and look whether the data are compatible with a straight line 
over a reasonable range of temperature. If this is the case the extrapolation to r^^^^ = 
should give an estimate for T^. Fig. ^ shows that for A = 2.0 the data is indeed rectified 
for 1.1 < T < 1.4, while outside of this temperature range the curves bend. Since plots 
for other reasonable choices of A look quite similar, the value of A can be estimated only 
within ±0.5. In all cases it is difficult to use the estimates for for finite N to extrapolate 
to the value of in the thermodynamic limit since the A^-dependence is rather weak and 
the error bars of Te,{N) are, due to the mentioned extrapolation, quite substantial. For the 
case of A = 2.0 a dependence of the form Td(A^) — Tc)(oo) oc 1/A^ seems, however, to be 
compatible with the data . 



In order to provide a more systematic way of extrapolating the relaxation times to the 



thermodynamic limit, we assume that the dynamical finite size scaling hypothesis |55| , [58|J59 
holds and make the Ansatz: 

T = N''f {N{T/Td - 1)^/^* } for A^ ^ oo and {T/Td - 1) ^ 0. (35) 

The scaling function f{Q must obey f{( oo) oc (~^* to recover the proper thermodynamic 
limit, i.e. Eq. (p4|). Using z* as a fit parameter (A is fixed to 2.0), we thus can try to generate 
a master curve from the t(T) curves for the different system sizes A^. That this is indeed 
possible if one chooses z* = 1.5 is shown in Fig. ^ From this figure we see that for large 
arguments the master curve does indeed show the expected power-law with an exponent —z* 
(dashed line). For T = Td the argument of f vanishes and hence we expect a A^-dependence 
of T of the form r oc A^^* and the inset of Fig. ^ shows that this is indeed the case. 

We mention that Eq. (|35D has a well-based theoretical foundation for second order phase 
transitions p2| , p5| , p8| , p9| , i.e. the case in which in Fig. ^ the temperatures Tq, Td, and 



Tg coincide at a unique critical temperature T^. The diverging relaxation time is then an 
immediate consequence of a diverging static susceptibility, and dynamic finite size scaling is 
a consequence of ordinary dynamic scaling [^. E.g., for second order transitions of mean- 
field spin glass models one has Eq. ( |55D with A/z* = •^mf + 2/3mf = 3, since the static 
mean field exponents of the spin glass order parameter and susceptibility are Pmf = 1 and 



•yMF = 1, respectively |2j]. Using the value of A = 2 ||22[, one hence finds for z* the value 
2/3. This result could be expected since the standard finite size scahng result for the critical 
relaxation in spin glasses with short range interactions is r oc with z = 4 in the mean 
field approximation ^2|]. (Here L is the linear dimension of the system.) This result can now 
be translated into the behavior of infinite range models at the marginal dimension d* = 6, 
i.e. at the dimension where mean field theory becomes valid, via A^ = L'^* , which yields 
T OC A^^/"^* = A^^*, i.e. z* = z/d* = 2/3. This result is also compatible with numerical 



simulations [BD|. However, the value z* ~ 1.5 found for the present model is clearly rather 



unusual and we are not aware of any analytical estimates for this exponent. 
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A further interesting question concerns the asymptotic decay of the correlation function 
C{t) towards qea as t — oo at T = Tp. In the context of the structural glass transition 
the time regime during which the correlation functions are close to the plateau is called the 

relaxation" whereas the decay below the plateau is called the "a-relaxation" 0. Mode- 
coupling theory predicts that at To the approach to the plateau is given by a a power law, 
i.e. 

Cit) - qEA oc r", T = Td. (36) 

A naive way to check for the presence of such a power law is to make a log-log plot of 
C{t) — qEA versus time. Fig. |10| shows such a plot for two relatively large systems (curves 
with open symbols) and we see immediately that there is no straight line, i.e. no power 
law dependence. However, one must recall that C{t) shows a significant dependence on N 
and that the prediction of Eq. ( |5B| ) seems to hold only for systems much larger than the 
ones studied here. Therefore we have tried to make an extrapolation, keeping time fixed, 
of the curves for finite N to determine the curve for N = oo. One possibility for such 
an extrapolation is to plot C{t,N) vs. While such a graph gives a straight line for 

> 480, for smaller A^ a curvature is clearly apparent. Therefore we tried whether a plot 
of C{t,N) vs. exp{— const ■ N), where the constant is a fit parameter, gives a straight line 
and found that for a broad range of times (from t = 10 to t = 300) this is indeed possible if 
the constant is approximately 1/400 

In Fig. |10| we have included the results of these two extrapolations also and we see that 
they do not give the same result. Since a priori it is not clear which type of extrapolation (if 
any!) is the correct one, it is difficult to tell what the shape of C{t) in the thermodynamic 
limit really is. It is very interesting to note, however, that the extrapolation with the 1/A^ 
dependence gives a curve C{t,N = oo) which is very well compatible with a power law of 
the form given by Eq. (^). Thus this gives some evidence that the extrapolation with 1/A^ 
is the correct one. The value for the exponent a we read off is 0.33 ± 0.04. 

It is also interesting to note that the theory predicts a one-to-one correspondence between 
the value of a and the exponent A from Eq. (|3^ ) [0. For a given value of a one can use 

r^(i + fe) ^ r^(i-^) 

r(l + 26) r(l-2a) ^ ^ 

to calculate b. (Here T{x) is the usual F-function.). The value of A is then given by 
l/2a + 1/26. If one uses the value A = 2.0 and the above relations one finds a = 0.36, in 
very good agreement with the value determined from Fig. |T0|. 

Two other important results from mode-coupling theory concern the shape of the correla- 
tors close to Td in the a-relaxation regime, i.e. in the time window where they fall under the 
plateau. The theory predicts that in this time regime the correlators can be approximated 
well with the Kohlrausch- Williams- Watts function, exp(— (t/r)^), a functional form that has 
been found to work very well in many glassy systems ||^,|33|j3^,|39| . We find, however, that 
even close to and the largest systems used, this functional form does not give a good fit 
to the data. 

The second prediction of the theory concerning the a-relaxation is the so-called time- 
temperature superposition principle. This principle implies that the correlator C{t) can be 
written as 
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C{t,T) = C{t/T{T)), (38) 

where C{x) is a temperature independent scaling function. The vahdity of Eq. (|38|) can be 
checked if one plots C{t,T) versus x = t/T(T). If the superposition principle is valid the 
curves for the different temperature should fall on a master curve for x ~ 1 and large x. 
For very small values of x, i.e. in the early /3-regime, no master curve is expected, since 



Eq. Pq ) is supposed to hold only in the a-regime. Fig. |ri| shows this kind of scaling plot 
and we see that even for a rather large system, N = 1280, there is no indication of such 
a time-temperature position principle. Of course, it is possible that this failure to verify 
Eq. ([38| ) is simply due to finite-size effects. Thus, it would be desirable to check Eq. (|38D 
for much larger systems. However, in view of the strong size effects on the relaxation time 



r near and below To, see Fig. ^ below, this is impossible for us with the present computer 
resources. 

Although we have just seen that the time-temperature superposition principle does not 
hold close to To for the accessible system sizes it is interesting to see whether this is the 
case also at lower temperatures. For N = 160 we have been able to go to temperatures as 
low as T = 0.7 and in Fig. |T2| we show the correlator as a function of t/r. If the curves 
for all temperatures are considered one finds that the superposition principle does again 
not hold (main figure). However, if one uses only the curves for the lowest temperatures, 
see inset, one finds that they all collapse onto a nice master curve. Thus we conclude that 
at sufficiently low temperatures the time-temperature superposition principle does indeed 
hold. We mention also, that the shape of this master curve is not an exponential, but that 
a stretched exponential with an exponent around 0.43 gives a satisfactory fit. 

In order to be able to understand this result a bit better, one needs to understand in more 
detail the relaxation of our model for T < To and finite A^, where all free energy barriers 
in phase space must obviously be still finite. One could argue that at low temperature the 
largest barrier dominates the dynamics and hence the relaxation depends on temperature 
only via a temperature dependent prefactor. This temperature dependence would have to 
be Arrhenius like and in order to check this we show in Fig. |T3] the T-dependence of r for 
the different system sizes. 

From this figure we see that at the lowest temperatures the T-dependence of r for 
the smallest system is indeed Arrhenius like. For temperatures around To and higher, 
one sees however significant deviation from this type of temperature dependence. Also 
for A^ = 320 one can see at the lowest temperature an Arrhenius law, but the activation 
energy is significantly larger than the one for A^ = 160. Since for increasing system size the 
lowest accessible temperature becomes higher and higher, it is not possible to see anymore 
the crossover from the non-Arrhenius T-dependence at intermediate temperatures to the 
Arrhenius dependence at small T. But the plot clearly shows that at To the relaxation 
times increases quickly with increasing system size, in agreement with the result from Fig. ^ 
(inset). Due to our present inability to equilibrate larger systems also significantly below To, 
we cannot determine reliably the A^-dependence of the activation energy of the Arrhenius 
law found at low temperatures. We mention that MacKenzie and Young found for the 
Sherrington-Kirkpatrick model , i.e. the mean-field Potts glass for p = 2, that for small 
systems (A^ < 128) and low temperature (T = 0.6Td = O.GTg) the relaxation times increase 



like t{N) oc exp{const ■ N^'^) |6^ whereas in a recent paper Billoire and Marinari |^ give 



evidence that the exponent of the power law is 1/3. If we consider a low temperature. 
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T = 0.7, this type of A^-dependence of r is compatible with our data with an exponent 0.5. 
However, if we determine the A^-dependence of the activation energy in the temperature 
regime where r(T, A^) shows an Arrhenius law, we find that this energy increases only very 
slowly, i.e. like log(A^) or a power of A^ with a small exponent (~ 0.1 = 1/p). (Note that 
the reason for the two different A^-dependencies is related to the fact that the prefactor of 
the Arrhenius law depends on A^ also.) 

Fig. O shows clearly that for small systems, A^ = 160, 320, and 640, it is possible to 
explore the region of temperatures below both the dynamical as well as the static transition. 
Note that in simulations of models for the structural glass transition it has never been 
possible to explore such low temperatures for comparable numbers of particles 0. On the 
other hand such models |^ do not seem to be much plagued by finite size effects, although 
for certain models for structural glasses they have recently been found |68 



V. RELAXATION OF INDIVIDUAL SPINS IN THE LOW TEMPERATURE 

PHASE 

In the previous section we have investigated the relaxation dynamics of the whole system 
and found that at low temperatures it shows a non-Debye behavior. In the present section 
we focus on the dynamics of the individual spins in order to obtain a better understanding 
for the occurrence of this non-exponentiality. 

In recent years it has been recognized that the non-exponential relaxation in supercooled 
liquids is often related to the so-called "dynamical heterogeneities" ||45| , |46[| . This means that 
the details of the relaxation dynamics of the individual particles (relaxation time, amplitude 
of the a-relaxation, etc.) is different for each different particle. The reason for this difference 
is (most likely) the fact that each particle has a slightly different neighborhood which thus 
affects the dynamics of the particle. Note that these differences are present only on the time 
scale of the a-relaxation r, since afterward the particle has changed its neighborhood and 
hence its characteristic dynamics. If the dynamics is averaged over a time much larger than 
r, all the particles behave the same. For spin glasses this is different, since the disorder is 
quenched. Hence the nature of the dynamics of the individual spins is an intrinsic property 
of each spin, since each spin is connected to the other ones by a set of different coupling 
constants. For a spin glass with short range interactions it is therefore not surprising that 
each individual spin has a different relaxation dynamics, and this is indeed what has been 
found in simulations For spins systems with long range interactions the presence of 



such dynamical heterogeneities is, however, not that clear, since each spin interacts with 
many different ones and hence one might argue that on average the different spins show the 
same relaxation dynamics. The goal of this section is to investigate this point in more detail. 

In order to characterize the dynamics of the individual spins we have calculated the 
autocorrelation functions Cj(t) for spin number i: 



^An exception are simulations of strong glass formers. E.g. in Ref. |64| it was shown that the 
system could be equilibrated even at temperatures as low as 0.8Td. However, one was still way 
above the Kauzmann temperature. 
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Q{t) = —-^{S,it')-S,{t' + t)). (39) 

Note that, in contrast to the case of structural glasses, it is here possible to average the 
right hand side over different time origins t', without loosing the information on the identity 
of the spin. Due to the single spin nature of the correlation function Ci{t) it is necessary to 
make this average over a sufficiently long time in order to obtain a reasonable statistics. We 
found that an average over 1000 a-relaxation times is needed, and therefore the following 
results have been obtained only for relatively small system sizes and 10 different samples for 
every temperature investigated. 

In Fig. [1^ we show the time dependence of Cj(t) for all the spins i = 1, . . . ,N for three 



different system sizes N. The temperature is Tq, i.e. the dynamical critical temperature 
at which the average relaxation dynamics, as measured by C{t), is already strongly non- 
exponential. From the figure we see that the relaxation dynamics for the different spins 
depends strongly on these spins in that, e.g., they relax to zero on time scales that span 
more than an order of magnitude. At a time where the correlation functions have reached 
0.5 of their initial value, the width of the range is even higher and increases rapidly with 
increasing system size. Furthermore we see from the figures that the curves for the individual 
spins seem to occur in clusters, i.e. that they do not fill the interval between the slowest 
and the fastest relaxation in a homogeneous way. Below we will discuss the reason for this 
clustering in more detail. 

In Fig. ^ we show the single spin autocorrelation function at a lower temperature. 
Comparing these curves with the ones in Figs. [1^ and |14|b we see that a decrease of T 



has made the distribution of the relaxation dynamics even wider. Also the presence of the 



clustering of the curves is now much more pronounced. From Fig. |T5| on also recognizes that 
the shape of the individual curves is not uniform at all since the ones which decay slowly 
tend to be, in the a-regime, much steeper than the ones that decay more rapidly. A more 
careful analysis shows that these slow spins show more or less an exponential relaxation 
whereas, as can already be seen from the figure, the fast ones show a strong deviation from 
a Debye law. Thus we conclude that the non-Debye behavior of C{t) found at low T, see 
Fig. is not due to a superposition of Debye laws with different relaxation times, but the 
sum of various different processes, some of which are Debye-like, some of which are not. 

In order to understand the microscopic reason for the presence of these dynamical het- 
erogeneities a bit better we have investigated to what extend the relaxation dynamics of an 
individual spin correlates with other quantities. For this it is necessary to characterize this 
dynamics in some way. As discussed above, the shape of the curves is not at all uniform, 
which makes such a characterization rather difficult. Therefore we decided to neglect all the 
variations of the shape completely and to characterize each curve just by the time it takes 
the spin to decay to a given value. Therefore we defined two different relaxation times, t^^'*^ 
and T^^''^\ via 

Q{t = rf ^^) = 0.4 and Ci{t = 7^^°■'^) = 0.7. (40) 



In Fig. 16 we show a scatter plot between (cj), the average energy of spin i, and the 
relaxation time, for both definitions of r^. We see that there is indeed a significant correlation 
between the energy and the relaxation time in that spins with high energy relax faster than 
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the ones with low energy. This result is very reasonable since a spin that has a low energy 
will be reluctant to change its value and therefore to go (with high probability) to a state 
with a higher energy. From the figure we also recognize that the correlation is present for 
both definitions of r^, from which we conclude that the details of this definition are not 
crucial. 



In order to investigate this point a bit closer we show in Fig. 17a a scatter plot of 



the relaxation time t-^'"^^ versus t-'^'^^ for the two temperatures. We see that although the 
correlation is not perfect, it is still very significant and therefore we conclude that the salient 
features of the correlation between the relaxation time and the mean energy shown in Fig. ^ 
will be observed even if a more careful characterization of the relaxation dynamics is made. 

Of further interest is the question how the relaxation time of a given spin at a given 
temperature is related to the relaxation time of the same spin at a different temperature. 
This dependence is related to the question of "chaos in temperature" , i.e. how the properties 
of a system change if temperature is changed. For mean field type system it is expected that 



these dependencies are rather weak [0,0 ■ agreement with this expectation we find that 
indeed the relaxation times for T = 0.9 are strongly correlated with those at T = 1.142, 
see Fig. p]7[b, irrespective of the definition of r^. Thus we see that this property seems not to 
be strongly affected by finite size effects. In passing we also mention that the mean energies 
(ej(T)) between the two temperatures are even much stronger correlated than the relaxation 
times [PH] . 

Before we end we come back to the observation discussed above that some of the single 
spin autocorrelation functions occur in clusters (see Fig. |15D. One potential reason that 
the relaxation dynamics of two spins is similar is that they are coupled together strongly, 
i.e. that their interaction constant Jij is large. In order to test this idea we identified for 
each realization of the disorder those spins that formed at T = 0.9 the cluster that relaxed 
slowest. (This identification was done visually by means of plots like the one shown in 
Fig. |TB|b). Say that this cluster involved k curves. We then determined the k{k — l)/2 
interaction constants between these k spins. The values of these constants are shown in 
Fig. 18 for ten different realizations of the disorder (filled circles). Also included in the 
figure is the Gaussian distribution of the coupling constants given by Eq. (^). From this 
figure we see that most of the points corresponding to the couplings are to the right of 
the mean of the distribution (vertical dashed line). Hence we conclude that the spins that 
form the slow cluster of relaxation curves are coupled together stronger than two arbitrary 
spins and therefore form a "dynamic entity". We note, however, that the fact that two 
spins are strongly coupled does not necessarily make them slow |51| which shows that such 
a strong coupling is only a necessary but not a sufficient condition for a slow dynamics. 

It is clear that the observations presented in this section are only modest first steps 
addressing the dynamics of the individual spins in the low temperature phase. It certainly 
would be interesting and useful to understand better how the distribution of the relaxation 
times of the spins depends on the system size and on temperature in order to obtain a better 
comprehension how the mean relaxation dynamics of the system is related to the one of the 
single spins. However, due to the high computational demand for this kind of investigations 
such studies have to be left to the future. 
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VI. CONCLUSIONS 



In this paper, we have presented the first detailed Monte Carlo investigation of the 10 
state mean field Potts glass model, for finite systems with sizes in the range from = 160 
to = 2560 spins. In the thermodynamic limit, it is known that this model exhibits both a 
dynamical transition at To where the system stops to be ergodic, and a static transition at 
To < Td where a glass order parameter q appears discontinuously (Fig. |^). The static spin 
glass susceptibility xsg remains finite both at Tq and T^. It would diverge only at a still 
lower "spinodal" temperature T, < Tq, if one were able to follow the disordered branch of the 
free energy at temperatures less than Tq. A further relevant temperature is the "Kauzmann" 
temperature Tk, Fig. |I]c, defined from the condition that the entropy of the high temperature 
phase vanishes. The spin autocorrelation function C [t) decays with time t at T > To in 
a two-step process, and the lifetime r of the "plateau" diverges as r oc (T — Td)~^, in the 
thermodynamic limit N —>■ oo (Fig. ^a). This behavior entails all qualitative features that 
one expects to be present at the structural glass transition (note that the equation of motion 
for C (t) indeed is described by the popular idealized mode-coupling theory proposed for the 
structural glass transition!). 

The questions asked in the present paper hence are as follows: can we verify these 
predictions from Monte Carlo simulations? How are the various transitions modified (i.e., 
rounded off) by the finite size of the considered model systems? Answers to these questions 
are not only of interest for a better understanding of the statistical mechanics of the present 
model system, they may also be useful to help with the interpretation of simulations of 
models for the structural glass transition. First of all, in the latter case the very existence of 
the various temperatures Td, Tq, Tg, Tk is still open to doubt. Secondly, even if one believes 
these temperatures should exist, their location for a particular model is still uncertain, unlike 
the present case where we have so much guidance from the exact solution. Of course, it is 
clear that a mean-field model is a rather special limit, and the sharpness of the dynamical 
transition at To probably is replaced by a smooth crossover from rather fast relaxation to 
very slow relaxation, as soon as one allows the interactions to become short ranged. In this 
sense, the finite mean-field Potts glass (where the singularity at To is rounded by the finite 
size of the system) may be qualitatively similar to the finite range model (although one 
should not push this analogy too far). 

Gratifyingly, we have established that the Monte Carlo results are qualitatively compat- 
ible with the theoretical predictions, although the finite size effects found were unexpectedly 
strong (i.e. they occur over a very wide temperature range as well) and they are not under- 
stood in detail, and hence we have found it too difficult to extract the various temperatures 
mentioned above directly from the simulation itself. E. g., for the sizes available, the method 
of looking for intersection points of the fourth oder cumulant or the Guerra parameter do 
not allow for a reliable estimation of Tq (Fig. Similarly, one could estimate the temper- 
atures Tg (Fig. H), Tk (Fig. |l]b) and To (Fig. ||) from a naive analysis of the data only very 
roughly. However, if one uses the theoretical knowledge on T^, one can estimate both the 
exponent A mentioned above and the exponent z* for the size dependence of the time r at 
Tc [t oc N^*^ from a dynamical finite size scaling analysis (Fig. We also found evidence 
that the predicted power law decay of the spin autocorrelation function C{t) — qea oc t~°' 
occurs for T ^ Tjj (Fig. |10]), but we could not confirm the expected time-temperature 
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superposition principle (Fig. |TT]). We also analysed the modified disconnected cumulants 
proposed in reference [|7ll, which should give a better estimation of the spin glass transition 
temperature in systems exibiting one step replica symmetry breaking patterns, but with 
the system sizes at our disposal they do not seem to work better than the corresponding 
connected parameters ||5T[] . 

Also some steps were taken to analyze the dynamics for T < To, by investigating the 
relaxation function Ci (t) of individual spins and corresponding relaxation times (Figs. 



171) . We find that the reason for the observed non-exponential relaxation of the mean 



relaxation function C{t) is related to the presence of a very strong dynamical heterogeneity. 
Furthermore we found that certain spins form dynamical clusters, the reason for which are 
likely their strong bonds between them. However, this mechanism seems not to be the only 
one and hence this point has to be investigated in the future in more detail. 

Thus, although many exact results are known for this model, and - unlike in models of the 
structural glass transition - we can equilibrate the system also at temperatures significanly 
below Td for a range of sizes (A^ < 640), we still are not able to answer many questions. 
Nevertheless, we think that the present model is a prototype model for glass transitions, 
and if better simulation algorithms become available, further studies of the present model 
should be very rewarding. 
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fig 1c 
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FIG. 1. (a) Energy e per spin plotted vs. inverse temperature (3 = 1/T for different system sizes 
(curves with symbols). The bold solid curve shows the one step replica symmetry breaking solution 
of DeSantis et al. [^3|,Q, the broken curve - which coincides with the former for T > Tq - is the 
replica symmetric solution, Eq. (25). The thin vertical lines indicated the inverse temperatures Pd 
(left) and Pq (right) of the dynamical transition and the static transition, respectively, (b) Analysis 
of the size dependence of the energy difference e]y{T) — eoo(T), using the one-step replica symmetric 
solution of DeSantis et al. |jl^ to calculate eoo(T). Inset shows the data of T = To = 1.142 plotted 
vs. A^"^/^ instead of N^^. (c) Entropy s per spin, normalized by its high temperature value, plotted 
vs. inverse temperature for different system sizes (curves with symbols). The bold dashed and 
the bold solid curve is the replica symmetric solution and the one-step replica symmetry breaking 
solution, respectivley. Vertical arrows indicate the static inverse transition temperature /3o and the 
inverse of the "Kauzmann temperature" Pk, where the entropy of the replica symmetric solution 
vanishes. fig2 




FIG. 2. Inverse of the reduced spin glass susceptibility XSG versus the square of inverse tem- 
perature for different system sizes (curves with symbols). The solid line shows the result of the 
replica-symmetric theory, Eq. (|2^). Inset: Plot of Xsg ~^ C^s/T)"^ to illustrate the nonmonotonic 
convergence towards Eq. (p6|). See main text for more details. 
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fig3c 




FIG. 3. (a) Order parameter distribution P{q) versus q for = 320 and various temperatures, 
(b) Order parameter distribution P{q) versus q for T = 0.7 and the three system sizes N = 160, 240 
and 320. The asymptotic value of the order parameter (from Ref. |5^) is included by a vertical 
line, (c) Value of the first moment J qP (q) dq of the order parameter distribution vs. A^~^/^. The 
inset shows that close to the transition temperature Tq ~ this moment scales like N~'^/'^. 
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fig4b 
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FIG. 4. (a) Fourth order cumulant (74 plotted versus temperature, for three values of N, 
N = 160, 320 and 640. The vertical straight line highlights the predicted static transition tem- 
perature Tq. (b) Same as (a) but for the Guerra parameter. The horizontal dashed line is the 
theoretical expectation for T <Tq. 
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FIG. 5. (a) Qualitative sketch of the mean-field predictions for the p-state Potts glass model 
with p > 4. The spin glass order parameter, qea, is nonzero only for T < Tq and jumps to zero 
discontinously at T = Tq. The spin glass susceptibility XSG follows a Curie- Weiss- type relation 
with an apparent divergence at < Tq, see Eq. (^). The relaxation time r diverges already at 
the dynamical transition temperature Tq. This divergence is due to the occurence of a long lived 
plateau of height qea in the time-dependent spin autocorrelation function C{t) . From Brangian 
et al. ||5^. (b) Temperature dependence of go and w(T) (see Eq. (|2^)) for p = 10, as obtained from 
the one-step replica symmetric solution of DeSantis et al. fl^ . 
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fig6b 
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FIG. 6. (a) Time dependence of the correlation function C{t), Eq. (13), for N = 1280 and 
various temperatures. Also included is data for the predicted values of the static, Tq, and dynamic, 
Td, transition temperature. The horizontal straight line shows the theoretical prediction from 
Ref. |l^ for the Edwards- Anderson order parameter at Tjj, qea = hm C{t), cf. Eq. p^). (b) Same 

— t^oo — 

as (a) but for the rotationally invariant correlation function Cmit) defined in Eq. (p^). 
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fig7b 
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FIG. 7. (a) Time dependence of the correlation function C{t) for T = 1.8 and for 
T = Td = 1.142 for several values of N. The solid line is the theoretical value of the Ed- 
wards-Anderson order parameter qEAiTn) for N ^ oo [13]. The dashed line shows the value we 

~7) . (b) Time dependence of the rota- 
Td = 1.142 for several values 
1.8. (c) 



use to define the relaxation time r. From Brangian et al. 
tionally invariant correlation function Cjii{t) for T = 1.8 and for T 
of N. The inset shows the limiting value CRi{t — > oo) as a function of N"^/"^ for T 
Parametric plot of CRi{t) vs. C{t) at T = To for different values of A^. The square indicates the 
plateau value obtained for A — > oo. The bold straight line describes the relation Cjij{t) = C{t), 
believed to hold for A ^ oo. 
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fig8 




FIG. 8. Temperature dependence of t~^^^ for different system sizes, using A = 2.0 as a trial 
value. The bold straight lines are fits on a proper subset of point. The resulting extrapolated 
values for T£){N) are quoted in the figure. 
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FIG. 9. Log-log plot of the scaled relaxation time t/N^* vs the scaled distance N{T 
from the dynamical transition temperature Tq, choosing z* = 1.5 and A/z* = 1.3. The inset is a 
log-log plot of r(r = Td) vs N. 
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figlO 




FIG. 10. Log- log plot of the relaxation function C{t) — qea versus time for T = Td, using 
the theoretical value of QeAj Eq. (^2]). The curves with the open symbols are the data from the 
simulation for two system sizes. The two curves with the filled symbols are the extrapolation of 
the simulation data to the case N = oo (see main text for details). The bold solid line is a fit with 
a power law. 

figii 




N = 1280. Temperatures from right to left:r = 1.8, 1.6, 1.5, 1.4, 1.360, 1.280, 1.240, 1.2, 1.142, 
and 1.131. The inset shows a magnification of the part of the curves for t/r > 1. 
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figi2 




FIG. 12. Plot of C{t) vs. t/r for N = 160. The temperatures are T = 0.7, 0.765, 0.85, 0.9, 
0.95, 1.0, 1.142, and 1.17 (left to right). The inset shows the same data, but including only the 
three lowest temperatures. 



fig13 




FIG. 13. Relaxation time r plotted vs. inverse temperature for different system sizes. The bro- 
ken vertical line indicates the location of the dynamical transition. Note the choice of a logarithmic 
scale for the ordinate. Error bars of r are mostly due to the sample-to-sample fluctuation. 
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fig14c 
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FIG. 14. Time dependence of the single spin autocorrelation function Ci(t) at T = To 
N = 160(a), N = 320 (b), and A'' = 640 (c). Each of the curves corresponds to a different spin. 
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fig15b 
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FIG. 15. Time dependence of the single spin autocorrelation function Ci{t) at T = 0.9 for 
N = 160 (a) and for = 320 (b). Each of the curves corresponds to a different spin. 
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FIG. 16. Scatter plot of Tj, the mean relaxation time of spin i as defined in Eq. (|40|), versus the 
mean energy (ej). The open and closed symbols correspond to T = 1.142 and T = 0.9, respectively. 
The points are for a typical sample of 320 spins. 
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FIG. 17. a) Plot of r^° '^^ [defined by Ci{t = n) = 0.7] versus rf ''^ [defined by Q(i = n) = 0.4], 

for T = 0.9 and T = 1.142 (open and filled symbols, respectively) showing that the two relxation 
times are correlated. Each point correspond to a different spin, b) Plot of the relaxation times 
rf and t^^''^'^ at T = 1.142 versus these relaxation times at T = 0.9. Each point correspond to 
a different spin. 
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fig 18 
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FIG. 18. Values of the bonds between spins with very slow relaxation in 10 different disorder 
realizations for 320 spins at a temperature T = 0.9 (filled circles). For clarity the points have been 
displaced vertically by various amounts. The continuous curve shows the Gaussian distribution 
from which the Jy are extracted and the vertical dashed line shows its mean. 
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